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In this article, we demonstrate that in a transport model of particles with kinetic con- 
straints, long-lived spatial structures are responsible for the blocking dynamics and the 
decrease of the current at strong driving field. Coexistence between mobile and blocked 
regions can be anticipated by a first-order transition in the large deviation function for 
the current. By a study of the system under confinement, we are able to study finite-size 
effects and extract a typical length between mobile regions. 
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out of equilibrium dynamics 

1. Large deviations in kinetically constrained models 
1.1. Large deviations formalism 

The large deviation theory can be viewed both as an extension and a new frame- 
work for developing statistical mechanics 1 , and relies on probabilistic foundations. 
This formalism yields information about the fluctuations of temporal trajectories 
in configuration space, in analogy with the usual canonical thermodynamics ap- 
proach which gives access to the fluctuations of observables at equilibrium such as 
the energy 

This approach can be extended to out-of-equilibrium systems, and in particular 
can be applied to Markovian dynamics where a well defined steady state exists: 
equivalents of free energies and entropies can be defined for observables extensive 
in time and can be computed either analytically (for very simple toy models) or nu- 
merically [2] . The tails of the distribution and the fluctuations of a given observable 
can then be quantified and reflect its sensitivity to the initial conditions. 



l 



March 3, 2013 16:23 WSPC/INSTRUCTION 

proceedings-PitardTurci'revision 



FILE 



2 F. Turd and E. Pitard 

In a very general way, for a given process having a well defined steady state, if 
O is a time-extensive observable we can write the associated partition function 

Z (s,t) = Yl Prob(history)e- s0(history) , (1) 

histories 

where a history is the time series of configurations {Ci, C2, . . . , Ct] and s is a 
transformation variable conjugated to the observable (via a Legendre transform), 
so that histories having a non typical value of O are selected according to the value 
of s. O can typically be the activity (the time integrated number of particle moves) 
or the integrated current, s can be viewed as a chaoticity temperature since it allows 
to select histories having either more and more rare events (for s positive) or more 
and more frequent events (for s negative). 

In the long time limit, the partition function behaves like 

Z {s,t) ~exp{t^o(s)}, (2) 

where ipo is called the large deviation function associated with the observable O. 

The function ipo(s) is the cumulant generating function of O, since all the 
cumulants of the underlying distribution can be obtained by means of derivation. 
For example the stationary value of O is simply given by 

Um i^Ml = -^ o (s)\ s=0 (3) 

t^-OO t 

where the average is taken over the ensemble of histories. Similarly, higher order 
cumulants are obtained from higher order derivatives. 

1.2. Application to kinetically constrained models 

In the context of the study of the glass transition and the jamming transition, the 
slow, frustrated dynamics of real glasses has been modeled under strong simpli- 
fications by particles/spin models where the moves/spin flips have some dynami- 
cal restrictions. Such toy models (referred to as kinetically constrained models or 
KCM [3j[4]) allow to explore elementary properties of glassy systems, such as long 
relaxation times, aging dynamics, anomalous diffusion etc. with a minimum number 
of ingredients. 

Given their relative simplicity, these models have been considered as good can- 
didates in the field of disordered systems for applying the concepts coming from the 
thermodynamics of histories. Many of the glassy properties related to the slowing 
down of the dynamics when the temperature decreases (or, equivalent ly, when the 
density of particles increases) are thought to be related to the formation of clusters 
of blocked particles slowly relaxing in the bulk, keeping the system amorphous on 
long time scales but collectively organized on shorter time scales: these patches of 
clustered particles, moving cooperatively in a non diffusive manner are often called 
dynamical heterogeneities [5}{7]. A natural time extensive observable allowing to 
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quantify the strength of such blocking and clustering processes is the total activity 
K(t), defined as the number of moves or spin flips from a reference time to = 
to time t: mobile patches will increase the total activity while blocked regions will 
keep it constant. 

The numerical study of KCM like the Fredrickson- Andersen model 8 , the Kob- 
Andersen 2d lattice gas etc. via the thermodynamics of histories can be made 



possible thanks to biased ensemble algorithms 10-12 : either in continuous or 
in discrete time, it is possible to make evolve a certain population N of different 
clones or replicas of the system and make them perform a biased dynamics so 
that well defined fluctuations of the activity are explored for a given value of s; 
such biased dynamics allows the direct computation of the cumulant generating 
function associated to the activity iJjk(s), which by construction is analogous to 
the Helmoltz free energy of equilibrium canonical ensembles. 

In the case of kinetically constrained models, the derivative of ^k(s) is discontin- 
uous at s = 0: this first order dynamical phase transition has been interpreted 13 
as the signature of separation in the space of possible histories between active histo- 
ries, where particles manage to escape from cages, and inactive ones, where particles 
are blocked by the dynamical restrictions. In particular, the first order transition 
at s = (unbiased dynamics) is interpreted as the coexistence of such histories in 
real space dynamics, this vision being compatible with the cage dynamics for which 
KCM have been proposed as model glass formers. 

What happens if the system is driven by an external force? Does the heteroge- 
neous behavior persist in terms of the large deviation approach? In order to address 
these questions, we consider a simple driven KCM model for granular materials. 



2. The Kob- Andersen model in presence of a drift 



M. Sellitto has recently proposed 14 an out of equilibrium variation of the original 
Kob- Andersen model, introducing a basic mechanism for mimicking an external 
shear over a granular fluid. This model can also be reinterpreted as the kinetically 
constrained version of a 2d Asymmetric Simple Exclusion Process (ASEP). 

sOn a 2d regular square lattice of size L x L, the Kob-Andersen model is a 
kinetically constrained model [9] where hard-core particles can move to any of the 
nearest-neighbour empty sites under the condition of having no more than two 
neighbouring particles before and after the move. The motion is isotropic, and the 
system is ergodic at any density, even if the relaxation times rapidly increase at 
high densities. 

Under periodic boundary conditions, it is possible to observe the set up of a 
macroscopic current as soon as an external bias is introduced. The simplest way to 
do it is to impose a uniform, unidirectional external field: following a discrete time 
dynamics, particles can hop in a given direction with probability p = min{l, exp(i5- 
df)} where df is the displacement vector from cell to cell. Combining these rules of 
motion and the Kob-Andersen kinetic constraint leads to the asymmetric exclusion 
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Fig. 1. (color online) a) Dynamic rules for the Kob- Andersen model in presence of a drift ; particles 
(red circles) can move and occupy empty sites if they have at least 2 empty neighbors before and 
after the move; an external field E acting on the whole volume is imposed in order to bias the 
dynamics in a given direction, b) Current vs field intensity relation for a system of linear size 
L = 100 for different values of the density; at low densities J(E) is a monotonic function, while 
at densities above p c ~ 0.79 c) the relation becomes non monotonic. 



processes with constrained dynamics [14] . 

At low densities the overall dynamics is well modeled by an effective 2d Asym- 
metric Simple Exclusion Process, without the Kob-Andersen constraints: the cur- 
rent is then proportional to the difference 1 — e _jE , which measures the difference 
between the hopping rates in the positive and negative field direction. This leads to 
a monotonic behavior for the current vs field relation J{E), which quickly saturates 
as the field intensity increases. We see that the introduction of the Kob-Andersen 
kinetic constraints does not lead to any important change in the dynamical pro- 
cesses with respect to normal ASEP (see fig.[l]b): this means that the correlations 
between particles movements do not go far beyond on-site hard core repulsion, at 
least for small densities. 

However, when the density of particles increases, the monotonicity breaks down. 
This happens approximately at density p c « 0.79 (fig. |^c) and the current eventu- 
ally saturates at a finite value at sufficiently large field intensity. The presence of 
a maximal value for the current allows to distinguish two characteristic dynamical 
regimes: for field intensities smaller than the optimal value E*(p) the relaxation 
times are comparable to what happens in the E = case, while for E > E* the 
dynamics slows down exponentially with the external field intensity [14] . 

Here we focus on the dynamics in the high density, strongly correlated regime, 
with particular attention to the sensitivity to initial conditions, the bimodality of 
the distributions of certain dynamical observables (in particular the activity and 
the current) and their scaling with the system size. 

2.1. Large deviations for the activity and the current 

Using the large deviation formalism introduced in section [l] both the activity (the 
number of microscopic movements in a whole history) and the total integrated 
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current have been studied 15 . From this analysis, it results 1) that a first-order 
dynamical transition exists for the activity for any value of the external field, being 
the signature of the presence of the kinetic constraints and 2) the current shows 
a first-order transition between an active and an inactive phase only in the limit 
of very strong fields. The interpretation of this phase transition occurring in the 
strong field regime (where the current decreases when the field increases) can be 
given by the accurate analysis of the changes that occur in configuration space when 
the external field is increased. 

The first order transition in the current large deviation function has to be in- 
terpreted as the coexistence of flowing and non flowing histories for the system, 
and a strong sensitivity to initial conditions. The flowing and non flowing histories 
phase coexistence is the signature of a bimodal distribution of the current in the 
thermodynamical limit. In practise (at s = 0), there is only an exponentially small 
number of trajectories that make the system blocked or with very small current; 
these rare histories are better sampled in the biased dynamics (at s > 0). We shall 
see in the following that this picture can be improved by a finite size analysis in 
the configurations space. 




Fig. 2. Numerical evaluation of the large deviation function for the current i/jq(s) for a system 
of linear size L = 60 at density p = 0.8 15 in the vicinity of s = 0. The discontinuity in the 
derivative of the large deviation function at the origin is clearly visible for strong fields. 



Indeed the origin of active and inactive patches resides on the microscopic dy- 
namics of the rearrangements of particles and the limitations due to the coupling 
of the kinetic constraint with the increase of the external field: these two ingredi- 
ents are responsible for a complete blockage of the dynamics in the limit of infinite 
lattice sizes and fields. Such microscopic rearrangements give rise to the transver- 
sal structures composed by filled and empty regions (that we have called domain 
walls) which allow for a corse-grained interpretation of the dynamics in terms of a 
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patchwork of active and inactive regions of quantifiable size. 



2.2. Blockage dynamics 

The variation of the external field intensity triggers modifications in the real space 
structure of the system. 

Once the external field is turned on, it is possible to observe 15 the formation 
of structures in the transverse direction: clusters of holes surrounded by bars of 
particles that eventually form T-junctions at the edges. These areas are not perma- 
nently blocked: defects (neighboring holes) are present so that the structures may 
relax, but the time needed for the relaxation is largely increased with respect to 
the undrifted case (fig. [3]). It is possible to measure the extension of the connected 
clusters of holes in the longitudinal and transversal direction. It appears (as dis- 
cussed in [15 ) that the transversal length is the most relevant and defines a walls' 
length. These walls of holes block the dynamics: their growth as a function of the 
external field is directly linked to the non monotonic behavior of the current as a 
function of the field. 
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Fig. 3. a)The square lattice Kob-Andersen models admits some frozen configurations composed 
by T-juctions of double- width occupied bars (particles in red). While these configurations are rare 
in the undriven case (E = 0), the application of an external field (E = 5 and p = 0.8 in panel 
b) ) enhances the probability of such configurations, where the residual mobility is assured by 
localized defects (holes) in the vicinity of the junctions, c) For a model of linear size L = 100 at 
density p = 0.8, connected clusters of holes allow the definition of a transversal and longitudinal 
length, wt and respectively. Transversal walls of holes (thick red line) grow when the external 
field increases while in the longitudinal direction (dotted line) weak field effects are measured. 
Average values have been computed over 100 different histories. 



3. Finite size analysis 

3.1. Bimodality of the current distribution in real space 

The role of transversal walls in the crossover from a monotonic to a non monotonic 
J(E) relation allows to shed light on the sensitivity to initial conditions and the 
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bimodal nature of the distributions of current and activity. 

We have chosen to modify the protocol of our simulations in order to explore the 
role of relevant lengthscales, in particular in the transversal direction. We still keep 
the torus geometry for our system (periodic boundary conditions in both direc- 
tions), but we now differentiate the transversal dimension H from the longitudinal 
dimension L: in order to keep the same statistics for particle moves for different 
system sizes, we rescale H in order to keep the product H x L ~ 60 2 constant, 
hence the number of particles is kept constant. 

Reducing the transversal size H makes the transversal blockages extension criti- 
cal for the current flow. If the blocked regions are of size comparable to the transver- 
sal dimension H, these regions slow down the flow in the system, leading it to an 
absorbing frozen state where no net flow is produced. The realization of such con- 
figurations depends both on the initial conditions and on the history. 




Fig. 4. Current vs field relation for different values of the thickness of the system H at density 
p = 0.80. Decreasing the transversal dimension (H < 30) increases the current fluctuations; 
although all averages are performed over 150 different histories, the noise of the data clearly 
increases when H decreases. 



Moreover, the growth of the fluctuations of the current is the signature of an 
important change in the shape of the distribution of the current. As shown in figure 
[5j large systems correspond to unimodal current distributions peaked around a well 
defined positive value. When the transversal size is reduced (to H = 20) a second 
peak appears in the vicinity of J = 0, whose relative amplitude increases when the 
system size is reduced. For very small H, this second peak at zero current becomes 
dominant, meaning that the histories corresponding to some non-zero flow become 
exponentially unlikely. 
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Fig. 5. Current distributions for strong fields (E = 2) at density p = 0.80 for various transversal 
dimensions H. Large systems have unimodal distributions peaked around a positive value of the 
current, while very small ones are peaked around zero. Intermediate sizes show bimodal distri- 
butions with different amplitudes of the peaks. The distributions are computed as kernel density 
functions over 100 different histories. 



3.2. Dynamical correlation length 

In order to better highlight the relation between structural properties and the 
dynamics, one can try to couple the changes in the dynamical properties as a 
function of the external field to the variation of the transversal system size H. To 
do so, we can use the definition of relaxation times given in [14] as the integral 

/>OC 

T rel {E,p,H)= / (/>(t)dt, (4) 



where <j)(t) is the persistence function quantifying the probability that the occupa- 
tion variable of a lattice site has never changed between time and t. As noted 
16 17 for the undrifted 2d Kob- Andersen model (which is ergodic for all densities 



below 1) it is possible to define a characteristic length scale 

-~ exp {i8(r^)}' (5) 

which gives an estimate of the spacing between mobile elements: in this sense it can 
be viewed as a dynamical correlation length. A system much larger than this length 
scale will have a well-defined relaxation dynamics. Below this length, the system 
may break into blocked patches that do not relax. In the case of a macroscopic flow, 
this breakdown of the system is clearly a limiting factor: no macroscopic current 
can pass through, reducing the current contribution to zero. Moreover, the smaller 
the system is, the longer will be the relaxation time, because more and more areas 
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will be frozen in their initial state. We compare the relaxation times of the system 
to its linear dimensions in order to extract the crossover length S p . In our case, this 
length will be dependent on the external field intensity E. 



p = 0.81 




10 20 30 

H 



Fig. 6. (colors online). Scatter plot of the relaxation times computed for different system sizes 
at different values of the external field. Each point corresponds to a unique history, so that for 
every couple (H, E) 20 different histories are considered. One can easily distinguish two groups of 
points: slow trajectories and fast ones, here distinguished by a gray dashed line. LOESS (LOcally 
Estimated Scatterplot Smoothing) best fits are also shown as dashed lines in order to illustrate 
the behavior of the average value of the relaxation times as a function of H. In the case of E = 
(black curve) the kink in the average corresponds to the value of the predicted critical length 
"0.81 ~ (orange dotted line), which rapidly grows for stronger fields and seems to saturate at 
a value S^ gl ~ 23(fuchsia line). 



The numerical study shown in figure [6] demonstrates how the sensitivity to initial 
conditions and the existence of a critical length are combined and depend on the 
intensity of the external field. In fact, while for sufficiently large systems one typical, 
almost constant relaxation time is well defined for any value of the external field, 
reducing the transversal dimension H leads to a rapid increase of the relaxation 
times at a certain H*(p,E); eventually very small systems are practically arrested. 
Moreover, around the critical length E) the distribution of the relaxation time 

splits in two parts, becoming bimodal: while a fraction of histories is characterized 
by the same relaxation time as large thermodynamic systems, a fraction of histories 
jumps to much higher relaxation times, getting trapped and slowed down by the 
formation of blocked patterns. 

If one takes the average relaxation times, one can observe that the kink cor- 
responding to the change of slope of the \og(r re i(H)} in the case of E — ap- 
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proximately corresponds to the critical length scale S p predicted in [17], so that 
we have H*(p, 0) ~ S p . The numerical simulations then suggest that when the 
field increases this critical length correspondingly increases: this is coherent with 
the picture provided before concerning the formation of new correlated structures, 
slowing down the dynamics and making harder and harder for particles to escape 
from their initial configuration. 

4. Conclusions 

We have discussed the dynamical implications due to an externally driven jamming 
transition for a simple 2d model. Dynamical heterogeneities have been related to 
the realization of rare rearrangements driven by the external field and the intrinsic 
bimodal nature of the distribution of current. A finite size effects analysis allowed 
to explore in real space such bimodality, previously signaled by a study in terms of 
large deviations functions. 
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